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We study how local equilibrium, and linear response predictions of transport coefficients are 
violated as systems move far from equilibrium. This is done by studying heat flow in classical 
lattice models with and without bulk transport behavior, in 1—3 dimensions. We see that linear 
response and local equilibrium assumptions break down at the same rate. The equation of state is 
also found to develop non-local corrections in the steady state. We quantify the breakdown through 
the analysis of both microscopic and macroscopic observables, which are found to display non-trivial 
f^*) ' size dependence. 

(N 

In studies of non-equilibrium systems, local equilibrium is an assumption which is essential to allowing the use of 
statistical mechanics and equilibrium or non-equilibrium thermodynamics Without local equilibrium, even the 
definition of temperature is not unique [Q, and it becomes unclear how to define simple transport processes. Local 
equilibrium is typically justified through conditions argued to be necessary or sufficient. This might involve checking 
1— 1 for Onsager reciprocity [gj, enforcing upper limits on local fluctuations of temperature (density,...) verifying that 
£^ I the equation of state holds locally ||, and so forth ||. However, these conditions and inequalities do not offer 
any quantitative guidance into how local equilibrium breaks down. Further, linear response has not been tested in 
conjunction with the breakdown of local equilibrium, which is an important consideration if one questions whether 
or not higher order corrections to Fourier's Law are consistent. Previous studies have observed the breaking of 
local equilibrium as well as deviations from linear response HJhJ. However few such cases are known, and the 
quantitative behavior of physical observables when local equilibrium is broken has not been studied previously. In 
||, the breakdown of local equilibrium was observed in the X-Y model and the Lorentz gas under thermal gradients, 
which was attributed to the infinite number of local conservation laws in the dynamics. Our results for non-integrable 
models will show that the integrability of the systems are not necessary for the deviations from local equilibrium to 
occur. In this letter we establish a quantitative guide to the rate at which concepts like local equilibrium and linear 
response become violated in systems which are subject to thermal gradients. We study lattice models in d = 1 — 3 
(/*■) ' spatial dimensions, including the Fermi-Pasta-Ulam (3 model, which does not have a bulk transport limit in d = 1. 
We also test how the equation of state is modified and its relation to expectations from irreversible thermodynamics. 

How should physical observables behave away from local equilibrium? A natural idea is that a physical observable A 
will deviate from its value in local equilibrium as we move further away from equilibrium. In our case, a temperature 
f| ■ gradient VT provides a natural measure of how far we are from equilibrium. Since the intrinsic physical properties 
should not depend on which side of the box is at a higher temperature, the deviation 5 A from its local equilibrium 
Ch 1 value is expected to behave as 



> 



h : 

While seemingly natural, such an analytic expansion is not trivial; in sheared fluids, transport coefficients have been 
seen to display non-analytic dependences on the shear rate, which have not been entirely clarified [jll). The coefficients 
C_4, C' A , . . . are in principle dependent on T and L, the size of the lattice in the direction of the gradient. If the relation 
is completely local, we expect that they will be independent of L. We shall find that the situation is more subtle. 
We study two systems which display qualitatively different transport behavior. Their Hamiltonians are 

H=\Y j [pI + ^^) 2 + v]. (2) 

r 

where V — (3(V(b r ) 4 /2 for the FPU— (3 model and V = </>r/2 for the <b 4 model. The (b 4 model has a well defined bulk 
limit for the thermal conductivity in d = 1 — 3 ficf| whereas the FPU model has L dependent thermal conductivity 
in d = 1 p^ , p^[ |. We will use (3 = 1 without loss of generality. Here r runs over all sites in the lattice (x = I, .... L, 
y,z = 1, JVj_), and the lattice derivative has components Vfc0 r = 4> r +e k ~~ 4>r ( e k is the unit lattice vector in the 
k-th direction). For d = 2, 3 we take iVi = 3 — 20 sites in the transverse directions with periodic boundary conditions. 
We thermostat the endpoints L = 0, N + 1 dynamically at temperatures T\ and T 2 , as discussed in O,[L0|. 
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Once we determine Ca for an observable A, it is possible to compute its non-equilibrium spatial distribution 
function. When n(T) behaves as a power law in the temperature range of interest, denoted n — cT^ 1 ', the temperature 
profile is T(x) = Ti(l - (1 - (T^/T^-^x/L) 1 ^ 1 -^ @. Such a power law behavior for k(T) has been shown to 
hold for the </> 4 model in d = 1 — 3 and also in most temperature regions in the d = 1 FPU f3 model, including the 
region we work with here Jl0|,[l3]]. The agreement for the predicted profile is shown in Fig. 0(a). Using Fourier's law 
(J = -«VT), Eq. (0) and T(x), we derive to leading order: 

^-Cj^)'-CJ T U' (3) 



.4. \ kT J \a + bx 

where a = 2\ 7 c/ J, 6 = 7 — 1, and J is the heat flow. So knowing C.4, we can also predict the spatial variation of 
the non-equilibrium observable A. This bring to light an interesting relation to coarse graining. Coarse graining in x 
over regions of length £, with L > £ > A, where A is the mean free path of the excitations, will provide no significant 
improvement towards recovering local equilibrium, since (up to the sign of CU), the functional form of ([?]) is positive 
definite. (This will be evident in Fig. 0(b).) 
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FIG. 1. (a) Temperature profiles for the (f) 4 theory and the FPU model and the corresponding normalized cumulants (b), for 
L = 162, with the boundary temperatures (0.1, 1.0). Symbols (o, x) are the analytic predictions using T(x) in (a) and Eq. (^) 
in (bh The different 7 between the two models in this temperature range causes a different relative sign between a and b in 
Eq. (W) accounting for the different shapes in (b). 



One might also be inclined to expand (|lj) in powers of (V"T)/T. However, in the region where the temperature 
dependence of the thermal conductivity can be described by a power law, k = cT~' 1 , one can show, using Fourier's 
law, that V™T/T — a n x (VT/T) n , where a n is a temperature independent constant. Strictly speaking, Fourier's law 
holds only close to local equilibrium, but as we shall see later, the deviations from it is of order (VT/T) 2 so that the 
difference is a higher order correction in the expansion in (VT/T). 

Local Equilibrium: With the local temperature given by T& = (p\), a natural measure for the deviations from local 
equilibrium is the deviation of the momentum distribution from the Maxwellian distribution. The cumulants {{p\)) = 

(Pk) - 3 (Pk) 2 > ((Pk)) = (Pk) " 15 (Pfc)(Pfc> + 30 (Pfc) 3 > and so on > normalized by the local temperature, {{p n ))/T n / 2 , 
provide a quantitative measure on how far we are from local equilibrium. In local equilibrium, ({p n )) = (n > 2). 
Consider first how systems typically behave under thermal gradients. In Fig. 0, the local temperature and ((p^,))/T 2 
are plotted against the position in both </> 4 and FPU systems. When the system is not too far from equilibrium, 
the temperature profile can be understood by locally applying Fourier's law J = -k(T)VT, where the heat flow J 
is constant |l^Jl^,|l^]. Since n depends on T, VT will also depend on T and hence the temperature profile becomes 
curved for increasing boundary temperature differences. We note that the temperature profile for the (j) 4 theory is 
visibly more curved under the same boundary conditions, even though the FPU is further from local equilibrium 
(Fig. |(b)h 

In Fig. [y(b), we see that the ((pf.)) cumulants are non-zero inside the system, so that we are no longer in local 
equilibrium. Contrary to naive intuition, the steepest gradient does not lead to the system being furthest from local 
equilibrium. In fact, the converse is true — the system is furthest from equilibrium in the flattest region. 

In Fig. ||, we plot the deviations of the 4-th cumulant of the momentum for the </> 4 theory and FPU model against 
VT/T, which corresponds to taking A = (p 4 ,) in ([j]). We are interested in the physics away from the boundaries 
and we shall always measure the physical quantities locally, well inside the system, although boundary effects can be 
readily understood |T^| . We find that VT/T provides a good measure of how far we are from equilibrium and in both 
models, the cumulants behave as 
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= c. 
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(4) 



where C LE = l.l(8)T°- 9 ( 2 ) (T = 1) and C[g u = 4.3(4)T°-"( 2 > (T = 8.8). Similar investigations at different T yields 
a weak T dependence for Cle which is difficult to establish. These results are consistent with d > 1 in both models at 
the same temperatures. Using C LE , C^ E U , we can predict the shape of ((p 4 ))/T 2 : In Fig. |(b), the non-equilibrium 
distribution ((0); symbols) is compared to simulation results (lines) agreeing nicely. 

Here, a relatively simple picture emerges: as we move away from equilibrium by increasing the difference in the 
boundary temperatures, each point in the interior deviates from local equilibrium in a predictable manner, without 
any threshold. Away from equilibrium, local equilibrium is an approximation that is quite good for small gradients 
since the deviations from it only vary as (VT) 2 . Similar results hold for higher momentum cumulants. 

The L dependence of Cj^ is quite intriguing. Naive argument suggests that since the gradients and the cumulants 
are local, the relation between them would not depend on L, at least in the (j) A theory where there is a bulk limit. This 
turns out not to be the case. In principle, it is possible that the effect we see will disappear in the large L, bulk limit. 
However, this seems implausible since we have excluded the region within the mean free path from the boundaries in 
the above results, using the properties of the model extracted in 10 

Linear Response: Let us now investigate the validity of linear response theory. This has been discussed previously as 
one of the criteria for the breakdown of local equilibrium even though no deviation was seen there. A priori, it is 
not clear if the linear response law can be used as a criterion for the breakdown of local equilibrium, rather than as 
an indication of higher order equilibrium corrections [T^ ] . Such an analysis assumes that even when linear response 
is broken, local equilibrium holds sufficiently well so that one can unambiguously define the temperature inside the 
system (or that we adopt a particular definition for T). The linear response prediction of the heat flow J is obtained 
by computing k from applying Fourier's law locally. This is denoted as Jlr, and agrees with direct measurements of 
J in the near equilibrium limit p0||. 
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FIG. 2. ({p 4 }) /T 2 away from equilibrium for the <fi 4 theory (d = 1, 2, 3, T = 1) and the FPU model (d = 1, 2) for various 
lattice sizes on an arbitrary scale. The fits of the data to a behavior const, x (VT/T) 2 are plotted for each lattice size. 
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FIG. 3. The deviations from the linear response law as a function of VT/T for the FPU model for T — 8.8 (left) and the 4 
theory for T = 1 (right) for d = 1. The dashes represent linear response and the solid line show the quadratic deviations from 
it of the form kVT [l + C LR {VT/T) 2 ] . 



When the temperature gra dient is small, linear response theory is applicable so that Fourier's law is satisfied globally: 
J = —n(T)(T 2 — T\)/L pc|Jl^| . As the gradient increases, curvature develops. When re ~ T~ 7 , we can integrate 
Fourier's law to obtain the next leading order correction due to curvature in T[x): 

Jlr - Jo = 7(7 + 1) a /VT X " 



This L 2 dependence simply indicates that Fourier's law is satisfied locally rather than globally. As the gradient 
increases even further, the energy that can be pumped through the system becomes less than that predicted by linear 
response theory even when it is applied locally ||lo|| . This is exactly the deviation we study here in Fig. ^|. This is 
quite difficult to measure, since unlike the cumulants, the theoretical value for the energy flow Jlr is not known and 
it needs to be obtained using linear response locally, which carries an error in itself. The relative deviation from the 
linear response result (as shown in Fig. |3|) can be reasonably well explained by 



5 LR = J - T ^ = C LR (^-\ , (6) 



^ 2 

Jlb. ~ V T 

where G^ R = -4(3)i 10 ( 2 ) (T = 1) and C[ PU = -6.6(8)L°' 9 W (T = 8.8). For large gradients, Eq. (g) will naively 
give rise to decreasing current with increasing gradient. However, when the gradient is this large, the higher order 
terms in (VT/T) becomes as important. We do not know how the system behaves under such an extremely large 
gradient, but it would be natural to expect that the current will saturate. The quadratic behavior in (||) is also seen 
for both the </> 4 theory and for the FPU model in d = 1 at other temperatures and also for the models in d > 1, even 
though the extraction of Clr involves larger errors in those cases. 

Within error, we see that the violation of linear response and local equilibrium are closely connected, occurring in 
the same manner: 

Sle ~ 5 LR . (7) 
Local equilibrium and linear response have no threshold, and break down at the same rate. 

Equation of State: The equation of state for these models arc simple in equilibrium since there is only one independent 
variable, which can be taken to be T. (In the FPU model in d = 1, there is also a possibly weak dependence on L in 
P eq ; unlike the L dependence of k, it is far less discernible.) We denote it P(T) (or E(T)) where P (E) is the pressure 
(energy density). We can measure them through the stress tensor, P — P xx — T 11 , E = 7~ 00 |k|. In Fig. ^, where 
we plot the relative deviation of the pressure from its equilibrium value, (P — P eq )/P eq , against VT/T. We see that 
the equation of state develops new dependences of the form 



P(T, VT, L) = P eq {T) 



1 + Cp 



(8) 



where C% = 1.5(1.2)i°- 9 ( 2 ) (T = 1), C£ PU = 4.1(6)L a30 ( 4 ) (T = 8.8). The non-equilibrium equation of state, 
P(T, VT, L) develops a non-trivial size dependence in Cp, rendering it non-local. Similar analysis for energy density 
yields the coefficients C% = 0.5(3)L 9 ( 2 ) (T = 1), C% PU = 1.7(7)L°- 3 M (T = 8.8). The quadratic behavior as in 
(^) is seen also in both models at different T for d = 1 — 3 but we do not have enough statistics to unambiguously 
extract Cp t E in those cases. While Extended Irreversible Thermodynamics (EIT) predicts the quadratic dependence 
in (||) for particle gases and liquids, precise identification of the non-equilibrium definitions of the physical quantities 
between the two theories is necessary before quantitative comparisons can be made. EIT predicts a local behavior for 
5P in contrast to our observations. 
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FIG. 4. Pressure P as a function of VT/T, away from equilibrium for ^ (top, T = 1) and FPU (bottom, T = 8.8). The 
fits of the data to a behavior const, x (VT/T) 2 are plotted for each lattice size. P is seen to increase away equilibrium in both 
models. (Vertical axis scaled as indicated.) 



We have quantified the violations of local equilibrium and linear response in </> 4 theory and the FPU (3 model in 
d = 1 — 3 dimensions, and observed that they break down at the same rate, with similar order of magnitude. Both 
are found to vary with the thermal gradient as (VT/T) 2 , and there is no threshold for the violations, appearing 
immediately as one moves away from global equilibrium. Other physical quantities such as the pressure and the 
energy density were also found to behave in a similar manner. We found that using the coefficients CU, we can 
predict the spatial dependence of non-equilibrium distributions of observables. As a consequence, coarse graining 
does not modify our conclusions. 

Since the definition of T is no longer unique when local equilibrium is broken, a question arises as to how the 
choice of non-equilibrium definition for T affects the results. Expressions for non-equilibrium deviations of J, P, E as 
in Eqs. (||),(§) will in general be affected covariantly; in particular, for a generic redefinition T = T' + v (VT/T) 2 , 
C' A = Ca + v{dA/dT). The local equilibrium violations seen in ((p™))/T™/ 2 as in Eq. (0) are invariant under such 
redefinitions, up to the order we consider. The physics, of course, is invariant under any redefinition in temperature. 
We find that the momentum cumulants provide the most natural and also the most clear criterion for the violation 
of local equilibrium. 

Certainly more questions remain: Importantly, we do not have an analytic understanding of how the coefficients 
Cle>Clr,Cp,Ce are related to the parameters of the theories, including their L dependence. Another question 
to consider is how other physical quantities behave away from local equilibrium. It would be interesting to further 
explore the consequences of non-locality in the equation of state, as well as use such models as a testing ground for 
non-equilibrium thermodynamics, since there is good control on the non-equilibrium steady state. 
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